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The results of subsonic jet validation cases for the Naval Research Lab’s Jet Engine 
Noise REduction (JENRE) code are reported. Two set points from the Tanna matrix, 
set point 3 (M. = 0.5, unheated) and set point 7 (M. = 0.9, unheated) are attempted 
on three different meshes. After a brief discussion of the JENRE code and the meshes 
constructed for this work, the turbulent statistics for the axial velocity are presented and 
compared to experimental data, with favorable results. Preliminary simulations for set 
point 23 (M. = 0.5, Tj/Tx = 1.764) on one of the meshes are also described. Finally, the 
proposed configuration for the farfield noise prediction with JENRE’s Ffowcs-Williams 
Hawking solver are detailed. 


Nomenclature 
L,Y, % Cartesian coordinates, with origin at the nozzle exit and centerline, 
x indicating the axial coordinate, increasing downstream 
r Radial coordinate, r? = y? + 2? 
A Approximate mesh spacing 
L Length of the mesh region meant to resolve the turbulent nearfield 
D; Exit diameter of the nozzle 
U; Fully-expanded jet exit velocity 
U, U’ Mean and root-mean square of axial velocity, respectively 
T; Fully-expanded jet exit temperature 
oe Mean and root-mean square of temperature, respectively 
To, T Startup and statistics-gathering simulation times, respectively 
0 Polar angle, with 6 = 0 pointing upstream (i.e. the negative-x direction) 
Yr Estimate of viscous length scale in the nozzle boundary layer 


I. Introduction 


Considerable effort has been invested by the computational aeroacoustic community in the prediction of 
subsonic jet noise through large eddy simulation (LES)#20 (see also the excellent review from Bodony24). 
Most of this work has been focused on using codes that employ high-order finite difference or finite volume 
schemes and structured grids in their discretization of the governing equations. A number examples of 
the application of unstructured LES codes to the problem of supersonic jet noise can also be found in 
the literature, e.g., the work done by researchers associated with Cascade Technologies22!241 and the Naval 
Research Laboratory (NRL) 2a83) Comparably little work has been done to validate the effectiveness of 
unstructured LES codes at simulating the acoustics of subsonic jets. While not needed to capture the 
geometry of the simple nozzles considered in this paper, the flexibility provided by unstructured meshes 
allows users to use small cell sizes in localized areas of the flow without affecting the mesh downstream, unlike 
structured meshes. Also, simulating the subsonic regime of jet flows approximating the takeoff condition is 
of interest to NASA’s Commercial Supersonic Technology (CST) project, with the goal of providing acoustic 
predictions through traditional Ffowcs- Williams Hawking techniqued24 and more sophisticated methods (e.g., 
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the forthcoming Generalized Acoustic Analogy work of Leib et al 85). The combination of these needs has 
motivated addressing this gap in the literature. 

The Jet Engine Noise REduction (JENRE) application developed by the NRL, is a computational 
aeroacoustics code designed to provide farfield acoustic predictions of jets using the LES approach to the 
simulation of turbulent flow. JENRE follows the Monotonically Integrated LES (MILES)241 philosophy of 
sub-grid scale modeling, relying on the inherent numerical dissipation of the code’s schemes to capture the 
effect of the unresolved turbulent scales will have on the resolved unsteady flow. Numerical dissipation 
is provided by blending low- and high-order finite element schemes through the Flux-Corrected Transport 
(FCT) approach. Proponents of this approach claim that the numerical dissipation effectively acts as an 
implicit model for the turbulent sub-grid scales, citing as evidence modified equation analysig38! and the 
results of applying such schemes to practical turbulent flows 222°] Numerous finite element schemes are 
available in the version of the JENRE code used in this work — here, the second-order Taylor-Galerkin 
schemd#4I#2] has been chosen for spatial and temporal discretization, coupled with a lumped-mass Galerkin 
scheme through the (FCT) method!##! 


II. Test Case Parameters 


The flow parameters for the simulations described in this work are taken from the Tanna matrix 4! and 
are shown in Table[I} The nozzle geometry chosen for this work was SMC000, the baseline round nozzle from 
the Small Metal Chevron series developed at the NASA Glenn Research Center (GRC). Described in more 
detail by Bridges and Wernet #4 SMC000 is an axisymmetric nozzle with a nominal two-inch exit diameter 
that has been subjected to extensive experimental investigation at the GRCBEF 


Set Point 3 | Set Point 7 | Set Point 23 
M 0.513 0.985 0.376 
My = U;/Ae0 0.500 0.900 0.500 
Ty Ts. 0.950 0.835 1.764 
po/D; 1.197 1.860 1.102 
Re; = 1" | 599-105 | 1.36-10° | 2.06- 108 


Table 1. Set points used in this work. 


Two mesh generation programs were used to create the meshes used in this work: Pointwis€ and 
Gmsh£2! The first two meshes, referred to here as PWO and PW1, were created with the Pointwise program, 
while the third mesh, GMO, was created with Gmsh. All meshes used the same SMC000 nozzle geometry. 
The boundaries of the nearfield, sponge zone, and outer domain were also identical. Figure |1] shows a 
diagram of the relative locations of the nozzle geometry, nearfield mesh region, and sponge zone and outer 
mesh boundaries. Diagrams of PW0, PW1, and GM0’s mesh spacing distributions are shown in Figures [2| 
and [4] respectively. Comparison of [2] and [3] shows that PW0 and PW1 differ in the element sizes in the 
region just downstream of the nozzle exit to about 12.5D;, with PW0’s spacing in the former axial location 
being twice as large as PW1’s, and increasing to the 0.04D; value more quickly. Figure [5] shows a crinkle 
slice of PW0’s z = 0 surface, where the surface of the slice has been constructed from the mesh element 
surfaces intersected by the z = 0 plane. Enlargements centered on the nearfield region for for this same 
crinkle slice are shown for the PWO and PW1 meshes, respectively, in Figures [9] and [10] The element edges 
are colored in blue in these two figures to make the differences in mesh density more apparent. 

Inspection of Figure[4]shows the GMO mesh used approximately the same mesh spacing sizes as the PW1 
mesh; however, the mesh spacing grows more gradually from the minimum 0.005D; to the maximum 0.04D; 
value. GMO also includes a refined region around the nearfield with target mesh spacing of about 0.05D,, 
intended to propagate acoustic disturbances to a Ffowcs-Williams Hawking integration surface. The impact 
of this change can be seen by comparing Figure [6] to Figure [11] 

The PWO and PW1 meshes contained approximately 15-10° and 27-10° nodes, respectively, while GMO, 
with its FW-H refinement region, has about 55- 10°. The nearfield region consists of a conical frustum of 
length L = 25D; with a diameter of D; at the nozzle exit and diverging at a rate of 4°. The shape of the 
nearfield was chosen to capture the development of the jet’s shear layer, and was informed by the shear layer 


2 of 22] 


American Institute of Aeronautics and Astronautics 


5D; <---> 


---->| 


sy al ae 


nearfield region 


25D; 


Figure 1. Diagram of the extent of the domain boundaries, including the size of the nearfield region and 
sponge zone. 
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Figure 2. Approximate mesh spacing distribution for the PWO mesh. 
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Figure 3. Approximate mesh spacing distribution for the PW1 mesh. 
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Figure 4. Approximate mesh spacing distribution for the GMO mesh. 
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Figure 5. z=0 slice of the PWO mesh. 
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Figure 6. z=0 slice of the GMO mesh, colored by instantaneous axial velocity. 


Figure 7. x = 2D; slice of the PW1 mesh, colored by instantaneous axial velocity. 
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Figure 8. x = 2D; slice of the PW1 mesh, colored by instantaneous axial velocity, enlarged figure. 
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Figure 9. Enlargement of the nearfield region of the z = 0 mesh slice shown in Figure [2] for the PW0O mesh, 
with cell edges highlighted to make the mesh spacing distribution more clear. 
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Figure 10. Enlargement of the nearfield region of the z = 0 mesh slice shown in Figure [3] for the PW1 mesh, 
with cell edges highlighted to make the mesh spacing distribution more clear. 


thickness growth rates reported by Bridges & Wernet 24 Figures [6| and seem to show that this choice 
was adequately large, as do Figures x = 2D; slices of the PW1 mesh. The nearfield length L is 
comparable tdL5IZ4] oy longer thar4BIM729/20| that used by other researchers, but the width of 25D; tan (4°) 
is considerably more narrow than similar simulation¢45/42129 found in the literature (4D, to 7.5D, is 
typical). 

For all meshes, anisotropic tetrahedra were extruded off the interior nozzle walls to partially resolve 
the boundary layer, likely to be quasi-turbulent according to Bridges and Wernet #8 For PWO and PWI1, 
this was done using Pointwise’s T-Rex tool, while Gmsh’s transfinite interpolation was used for GMO. The 
spacing of these cells in the wall-normal direction is initially 5-10~4D,, and grows at a rate of 1.2 until 
the tetrahedra become isotropic. The length scale of the viscous sublayer for the flows studied here is 
approximately y; = 4.1-10~°- D, (estimated from smooth flat plate boundary layer theory as described by 
Schlichting#%), and thus the simulations will not meet the “wall-resolved” requirement of A = 2y, 20 

The minimum mesh spacing in the nozzle nearfield region is 0.01D,; for the PWO mesh, and 0.005D; for 
the PW1 and GMO. Other subsonic jet LES examples in the literature use minimum mesh spacings ranging 
from about 0.002D BEF to 0.003.D, E54 and past supersonic jet calculations done by NRL researchers used 
spacings as small as 0.0035D,; in shear layer2%24 The mesh spacing increases gradually downstream to 
0.04D; at the end of the nearfield region for both meshes. 

Boundary conditions consisted of adiabatic no-slip for all nozzle walls, and the appropriate constant 
stagnation pressure and temperature for the set point under investigation for the nozzle inflow. JENRE’s 
characteristic-based farfield boundary condition was used for all other boundaries. As indicated in Figure[I] 
a sponge zone was used to damp any spurious reflections from these farfield. 


III. Results 


A. Set point 3 


For set point 3 and the two Pointwise meshes, a constant time step was chosen such that the local CFL number 
would not exceed 0.65 anywhere during the calculation. This time step size was reevaluated periodically 
during the run (about every 100 steps). The GMO runs used a constant timestep of 5.13 - 10~°D,;/Uj, 
corresponding to a CFL of about 0.24. Each set point 3 run was started from a uniform flow of zero velocity 
and ambient pressure and temperature. Table |2|shows the amount of time the simulation was run before 
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| To -U;/D; | T-U;/D; 


PWO 215.0 300.0 
PW1 747.0 547.0 
GMO 202.0 145.0 


Table 2. Set point 3 startup (70) and statistics (T) times. 


gathering statistics (Jy) and during the statistics-gathering process (J). The values for T - U;/D, for the 
PWO and PW1 meshes comparable to the periods of 87.5D;/U; 2% 500D,; /U; = and 833D,;/U/ found in 
the literature. GMO, the largest and thus most computationally expensive mesh, may benefit from more 
time. 


Figure 11. Instantaneous axial velocity non-dimensionalized by U; for set point 3 run on the PW1 mesh. 


After performing the simulations with the JENRE code as described in the previous section, the resulting 
data was post-processed to obtain mean and RMS statistics of axial velocity, and then compared to Bridges 
and Wernet’s measurements! Figure[11]shows the instantaneous axial velocity from an arbitrary time step 
for set point 3 on the PW1 mesh, and seems to indicate that the width of the mesh’s nearfield region is large 
enough to capture the majority of the turbulent flow field. 

Figure[12] and [13|]show the mean and RMS of axial velocity along the centerline of the jet as a function of 
the distance from the nozzle exit, and compare these values to the “consensus” dataset found in reference #6 
which consists of a representation of an expected value for these quantities informed by experiments run by 
many researchers performed at multiple labs over a considerable length of time, including to a significant 
degree by Bridges and Wernet themselves. 

Figure shows that little difference is observed between the results for the mean axial velocity U 
along the centerline for the three meshes. All simulations predict the length of the potential core (the region 
downstream of the nozzle exit that is unaffected by the shear layer) accurately, and consistently underpredict 
U downstream slightly. More spread between the two simulations is seen in Figure [13}s root-mean square 
axial velocity U’ results. Both simulations match the experimental data reasonably well for « < 10.D;, and 
show significant dissipation for the region beyond about z > 15D,. Interestingly, the simulation performed 
on the coarser PWO mesh agrees better than the the two finer meshes. GM0’s RMS results show a bit more 
variation over x > 10Dj;, likely a result of the shorter statistics sampling time for this run. 
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Figure 12. Mean axial velocity U along centerline (r = 0) for set point 3. 
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Figure 13. Root-mean square of axial velocity U’ along centerline (r = 0) for set point 3. 
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Figure 14. U along lipline (r = D;/2) for set point 3. 
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Figure 15. U’ along lipline (r = D;/2) for set point 3. 
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The U and U’ data along the lipline (i.e., r = 0.5D,) as a function of axial distance from the nozzle are 
shown in Figures[14] ana [15] respectively. Like the U centerline results in Figure [12] there is little difference 
between the three simulations, and all simulations tend to underpredict Bridges & Wernet’s measurements, 
with the largest discrepancy just downstream of the nozzle exit. Figure shows the simulations also 
have trouble in the same region for the U’ data, where a noticeable overshoot in the JENRE predictions is 
observed. This spike in U’ is often observed when the boundary layer unsteadiness in the simulations is too 
low, leading to the periodic shedding of a large, laminar-like vortex. The magnitude of the spike observed 
here is smaller than that of other simulations, however. 


B. Set point 7 


The set point 7 runs, like those of set point 3, were started from uniform flow of zero velocity and ambient 
pressure and temperature, with the exception of the run on the GMO mesh, which was started from the 
end of the corresponding set point 3 run in an attempt to decrease the amount of time needed to reach a 
statically steady state.. Table[2]shows the amount of time the simulation was run before gathering statistics 
(To) and during the statistics-gathering process (JT). With the interest of decreasing the turnaround time of 
the simulations, the Ty values are lower than those used for the set point 3 runs. 


| To -U;/D; | T-U;/D; 


PWO 69 170 
PW1 67 182 
GMO 45* 63 


Table 3. Set point 7 startup (To) and statistics (T) times. The GMO run was started from the set point 3 flow 
at 166D;/U;. 


Like the set point 3 simulations, the PWO and PW1 runs used a constant time step that prevented the 
CFL number from exceeding 0.5 for the startup period. This approach was also used for the GMO run, 
but a lower CFL limit of 0.25 was required to maintain stability. For the statistics-gathering period T, a 
constant time step was used for all simulations: 1.23 -10~4D,;/U; for PW0, 1.29-10~*.D,;/U; for PW1, and 
6.95 -10~-°D,;/U; for GMO. 

Figures and show the axial velocity mean and RMS results, respectively, for the centerline. The 
results are very similar to the corresponding set point 3 plots. The potential core is well predicted by all 
three simulations, and the mean flow after the potential core is slightly under-predicted, with satisfactory 
agreement overall. The location of the peak RMS is a bit upstream of the experimental data, but not as 
much as the set point 3 results. 

Figures [18] and [19| show the axial velocity mean and RMS results, respectively, for the r = 0.5.D, lipline. 
Like the centerline results, the data is very similar to set point 3. The mean axial velocity is underpredicted 
just downstream of the nozzle, but corresponds well to the experimental data otherwise. There is a small 
RMS spike just downstream of the nozzle exit, and the RMS data is lower than the experiment for x > 12D,;, 
and may benefit from more simulation time, especially the GMO mesh. 


C. Set point 23 


The preliminary set point 23 run was initiated from the last time step of the set point 3 calculation, with 
the goal of reducing the amount of time needed for the flow to reach a statistically stationary state. This 
set point has only been performed on the PW1 mesh to date (the finest mesh available at the time the 
simulations were started). The JENRE data was post-processed for the same U and U’ statistics presented 
in the setpoint 3 section, and also for the mean and RMS temperature data (T' and T’, respectively) along the 
nozzle centerline and lipline. The temperature statistics were compared to Rayleigh scattering measurements 
performed by Mielke & Elam 21 

Figures and show U and U’ for the setpoint 23 case, respectively. The agreement between the 
JENRE and experimental data is striking in the U plot. The U’ results in Figure [21|do not show the same 
excellent correspondence with the measurements, but still agree fairly satisfactorily with the experimental 
data. Interestingly, for the c > 10D, region, JENRE’s U’ prediction is better than what is observed for 
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Figure 16. Mean axial velocity U along centerline (r = 0) for set point 7. 
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Figure 17. Root-mean square of axial velocity U’ 
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Figure 18. U along lipline (r = D;/2) for set point 7. 
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Figure 19. U’ along lipline (r = D;/2) for set point 7. 
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Figure 20. U along centerline (r = 0) for setpoint 23. 
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along centerline (r = 0) for setpoint 23. 
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setpoint 3 in Figure[13| Comparison between the U’ experimental data for the two set points seems to show 
that the heated jet case tends to have lower axial velocity fluctuations in the x > 10D; region, perhaps due 
to the more intense mixing associated with the hot jet’s temperature gradients. 
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Figure 22. U along lipline (r = D;/2) for setpoint 23. 


The U and U’ lipline data for the setpoint 23 simulation are shown in Figures and [23} respectively. 
Again, JENRE’s prediction for U closely approximates Bridges and Wernet’s measurements, and the U’ 
agreement is also fairly good. Interestingly, the overshoot just downstream of the nozzle exit seen in setpoint 
3’s U’ lipline data Figure[15) is not observed for this flow condition. 

eee mean and RMS temperature data (T and T’, respectively) along the centerline, 
along with the Rayleigh scattering measurements from Mielke and Elam2! The T data agrees well for 
x < 7D,;, but diverges from the measured values for x > 7D,;, where it is a bit higher than the experiment. 
Similarly, the T’ data in Figure [25] agrees well for « < 10D,, but not for the points further downstream. 
(Mielke cautiong=2! that the data downstream of x = 10D; my be contaminated by particles from the ambient 
flow, which would interfere with the Rayleigh scattering technique used for the temperature measurement. 

The T and T” statistics for the lipline region for the present setpoint are shown in Figures and 
respectively. JENRE’s T data follows the experimental trend fairly well, but is shifted up for the entire 
axial range. Finally, the T’ data is somewhat close to Mielke & Elam’s measurements in the region just 
downstream of the nozzle exit, but does not compare well as x increases. Again, Mielkd®=2! indicates a higher 
uncertainty for the data in this region. More investigation is needed. 


IV. Conclusions and Future Work 


This paper has presented the preliminary results of applying the JENRE code to the task of predicting 
the nearfield turbulence of three subsonic jets. Three meshes of the SMC000 geometry were generated and 
then used with the JENRE LES code to simulate set points 3 and 7 from the Tanna test matrix, with 
one mesh also used to simulate set point 23. Good agreement between the JENRE data and experimental 
measurements was observed for the axial velocity statistics in the nearfield, indicating that JENRE holds 
promise as a LES/CAA tool. JENRE’s results for the mean axial velocity along both the lipline and centerline 
closely followed the experimental data for all three set points. The RMS axial velocity agreement was also 
quite good, but may benefit from a bit more simulation time. Interestingly, little difference between the 
results for PW0O (the coarsest mesh, with a minimum shear layer mesh spacing of 0.01D;) and PW1 and 
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Figure 23. U’ along lipline (r = D;/2) for setpoint 23. 
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Mean temperature T along centerline (r = 0) for setpoint 23. 
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Figure 25. RMS of temperature T’ along centerline (r = 0) for setpoint 23. 
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Figure 26. T along lipline (r = D;/2) for setpoint 23. 
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Figure 27. T’ along lipline (r = D;/2) for setpoint 23. 


GMO (minimum shear layer spacing of 0.005.D;) was observed for the turbulent statistics presented here — 
though it remains to be seen if this is true for the planned acoustic predictions. 

In the future, JENRE’s Ffowcs-Williams Hawking (FW-H) solver will be used to obtain farfield acoustic 
predictions. Figure shows the proposed location of the FW-H integration surface superimposed on the 
non-dimensional vorticity magnitude for the PWO mesh. Following the advice of Erwin and Sinha23) the 
surface has been placed where the vorticity falls below 0.5, which should hopefully keep the integration 
surface outside of the region where hydrodynamic pressure fluctuations are present. For this first attempt 
at farfield noise predictions, no endcaps or disk averaging will be used at the FW-H surface outflow — it is 
hoped that the surface will be long enough to contain all the significant jet noise sources. 
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